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Abstract. Since 1997 the afterglow of 7-ray bursting sources (GRBs) has occasion- 
ally been detected in the radio, as well in other wavelengths bands. In particular, the 
interesting and unusual 7-ray burst GRB980425, thought to be related to the radio 
supernova SN1998bw, is a possible link between the two classes of objects. Analyzing 
the extensive radio emission data available for SN1998bw, one can describe its time 
evolution within the well established framework available for the analysis of radio emis- 
sion from supernovae. This then allows relatively detailed description of a number of 
physical properties of the object. The radio emission can best be explained as the inter- 
action of a mildly relativistic (r ~ 1.6) shock with a dense preexplosion stellar wind- 
established circumstellar medium (CSM) that is highly structured both azimuthally, 
in clumps or filaments, and radially, with observed density enhancements. Because of 
its unusual characteristics for a Type Ib/c supernova, the relation of SN1998bw to 
GRB980425 is strengthened and suggests that at least some classes of GRBs origi- 
nate in supernova (SN) massive star explosions. Thus, employing the formalism for 
describing the radio emission from supernovae (SNe) and following the link through 
SN1998bw/GRB980425, it is possible to model the gross properties of the radio and 
optical/infrared (OIR) emission from the half-dozen GRBs with extensive radio obser- 
vations. From this we conclude that at least some members of the "slow-soft" class of 
GRBs can be attributed to the explosion of a massive star in a dense, highly structured 
CSM that was presumably established by the preexplosion stellar system. 



1 SN1998bw/GRB980425 

The suggestion of an association of the Type Ib/c SN 1998bw with the 7-ray burst 
GRB980425 may provide evidence for another new phenomenon generated by 
SNe - at least some types of GRBs may originate in some types of SN explosions. 
Because SN1998bw/GRB980425 is by far the nearest and best studied of the 
7-ray bursters, it is worthwhile to examine its radio emission in detail before 
proceeding to the discussion of the radio emission from other GRBs. 



1.1 Background 

While generally accepted that most GRBs are ex trem ely distant and energetic 
(see, e.g., the discovery of GRB980425 |07) on 1998 April 25.90915 

and its possible association with a bright supernova, SN1998bw at RA(J2000) 
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= 19 h 35 m 03?31, Dec(J2000) = -52°50'44"7 pH , in the relatively nearby spi- 
ral galaxy ESO 184-G82 at z = 0.0085 (distance ~ 40 Mpc for H = 65 



km s" 1 Mpc" 1 ) p|gM|6p,114,128|, has introduced the possibility of a 



SN origin for at least some types of GRBs. The estimated explosion date of 
SN1998bw in 1998 between April 21 - 27 @ corresponds rather well with the 
time of the GRB980425 outburst. Iwamoto et al.[[32| felt that they could restrict 
the core collapse date for SN1998bw even more from hydrodynamical modeling 
of exploding C + O stars and, assuming that the SN1998bw optical light curve 
is energized by 56 Ni decay as in Type la SNe, they then placed the coincidence 
between the core collapse of SN1998bw to within +0.7/— 2 days of the outburst 
detection of GRB980425. 

Classified initially as an SN optical Type lb [pif , then Type Ic then 
peculiar Type Ic [E3lM, then later, at an age of 300 - 400 days, again as a 



Type lb 83 , SN1998bw presents a number of optical spectral peculiarities that 
strengthen the suspicion that it may be the counterpart of the 7-ray burst. 

When the more precise Beppo-Sax NFI was pointed at the Beppo-Sax error 
box 10 hours after the detection of GRB980425, two X-ray sources were present 
[ jS4| . One of these, named SI by Pian et al.Q), was coincident with the position 
of SN 1998bw and declined slowly between 1998 April and 1998 November. The 
second X-ray source, S2, that was ~ 4' from the position of SN 1998bw, was 
not (or at best only marginally with a < 3er possible detection six days after the 
initial detection) detectable again in follow up observations in 1998 April, May, 
and November Q|5) . 

However, although concern remains that the Pian et al.j84|,^5| X-ray source, 
S2, might have been the brief afterglow from GRB980425 rather than the Pian 



source SI associated with SN1998bw, Pian et al. 1 85 concluded that SI has a 
"high probability" of being associated with GRB980425 and that S2 is more 
likely a variable held source. 



1.2 Radio Emission 

The radio emission from SN 1998bw reached an unusually high 6 cm spectral 
luminosity at peak of ~ 6.7 x 10 28 erg s _1 Hz -1 , i.e., ~ 3 times higher than 
either of the well studied, very radio luminous Type Iln SNe SN 1986J and SN 
1988Z, and ~ 40 times higher than the average peak 6 cm spectral luminosity 
of Type Ib/c SNe. It also reached this 6 cm peak rather quickly, only ~ 13 days 
after explosion. 

SN 1998bw is unusual in its radio emission, but not extreme. For example, 
the time from explosion to peak 6 cm luminosity for both SN 1987A and SN 
1983N was shorter and, in spite of the fact that SN1998bw has been called "the 
most luminous radio supernova ever observed," its 6 cm spectral luminosity at 
peak is exceeded by that of SN 1982aa |129[ |. However, SN 1998bw is certainly 
the most radio luminous Type Ib/c radio supernova (RSN) observed so far by a 
factor of ~ 25 and it reached this higher radio luminosity very early. 
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1.3 Expansion Velocity 

Although unique in neither the speed of radio light curve evolution nor in peak 
6 cm radio luminosity, SN1998bw is certainly unusual in the combination of 
these two factors - very radio luminous very soon after explosion. Kulkarni et 
al.|]6^| have used these observed qualities, together with the lack of interstellar 
scintillation (ISS) at early times, brightness temperature estimates, and physical 
arguments to conclude that the blastwave from SN1998bw that gives rise to the 
radio emission must have been expanding relativistically. On the other hand, 



Waxman & Loeb 1 121 argued that a sub-relativistic blastwave can generate the 
observed radio emission. However, both sets of authors agree that a very high 
expansion velocity (;> 0.3c) is required for the radio emitting region under a 
spherical geometry. 

Simple arguments confirm this high velocity. To avoid the well known Comp- 
ton Catastrophe, Kellermann & Pauliny-Toth |5^] have shown that the bright- 
ness temperature Tb < 10 12 K must hold and Readhead |88| has better defined 
this limit to Tb < 10 11 ' 5 K. From geometrical arguments, such a limit requires 
the radiosphere of SN1998bw to have expanded at an apparent speed ;> 230, 000 
km s _1 , at least during the first few days after explosion. Although such a value 
is only mildly relativistic (T ~ 1.6; T = (1 — ^)~*), it is still unusually high. 
However, measurements by Gaensler et al.[p4| and Manchester et al.|Q have 
demonstrated that the radio emitting regions of the Type II SN1987A expanded 
at an average speed of ~ 35,000 km s" 1 (~ 0.1c) over the 3 years from 1987 
February to mid-1990 so that, in a very low density environment, such as one 
finds around Type Ib/c SNe, very high blastwave velocities appear to be possible. 



1.4 Radio Light Curves 



An obvious comparison of SN1998bw with other radio supernovae (RSNe) is the 
evolution of its radio flux density at multiple frequencies and its description by 



known RSN models (see [ 126 1 ) . The radio data available at tittp://www.narrabri. 
at nf . csiro . au / ~mwier ing /grb/grb980425/ are plotted in Fig. |1|. SN1998bw shows 
an early peak that reaches a maximum as early as day 10 - 12 at 8.64 GHz, 
a minimum almost simultaneously for the higher frequencies (y > 2.5 GHz) 
at day ~ 20 - 24, then a secondary, somewhat lower peak after the first dip. 
An interesting characteristic of this "double humped" structure is that it dies 
out at lower frequencies and is relatively inconspicuous in the 1.38 GHz radio 
measurements (see Fig. |l]). 

Such a "double humped" structure of the radio light curves can be reproduced 
by a single energy blastwave encountering differing CSM density regimes as 
it travels rapidly outward. This is a reasonable assumption because complex 
density structure in the CSM surrounding SNe, giving rise to structure in the 
radio light curves, is very well known in such objects as SN1979C [76 122,123], 
SN1980K f Ugl, a nd, particularly, SN1987A g|. 

Weiler et al. 126 pointed out what has not been previously recognized, that 
there is a sharp drop in the radio emission near day ~ 75 and a single measure- 
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Fig. 1. The radio light curves of SN 1998bw at 3.5 cm (8.6 GHz; upper left, open 
circles, solid line), 6.3 cm (4.8 GHz; upper right, stars, dashed line), 12 cm (2.5 GHz; 
lower left, open squares, dash-dot line) and 21 cm (1.4 GHz; lower right, open triangles, 
dotted l ine) . The curves are derived from a best fit model described by the Equations 1 
- 10 of [126] and the parameters and assumptions listed in Table |l|. During the 50 day 
intervals from da y 25 - 75 and from day 165 - 215 the emission and absorption terms 
(Ki and K3; see [126]) increase by factors of 1.6 and 2.0, respectively, corresponding 



to a density increase of 40% with a 6 day boxcar smoothed turn-on and turn-off of the 
enhanced emission/absorption. 



ment epoch at day 192 that is significantly (~ 60%) higher at all frequencies 
than expected from the preceding data on day 149 and the following data on 
day 249. 

Weiler et al. were able to explain both of these temporary increases in radio 
emission by the SN blastwave encountering physically similar shells of enhanced 
density. The first enhancement or "bump" after the initial outburst peak is 
estimated to start on day 25 and end on day 75, i.e., having a duration of 
~ 50 days and turn-on and turn-off times of about 12 days, where the radio 
emission (Ki) increased by a factor of 1.6 and absorption (K 3 ) increased by a 
factor of 2.0 implying a density enhancement of ~ 40% for no change in clump 
size. Exactly the same density enhancement factor and length of enhancement is 
compatible with the "bump" observed in the radio emission at day 192 (i.e., the 
single measurement within the 100 day gap between measurements on day 149 
and day 249), even though the logarithmic time scale of Fig. [l] makes the time 
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interval look much shorter. The decreased sampling interval has only one set 
of measurements altered by the proposed day 192 enhancement, so Weiler et 
al.[126| could not determine its length more precisely than < 100 days. 

Li & Chevalier |35| proposed an initially synchrotron self-absorbed, rapidly 
expanding blastwave in a p oc r~ 2 circumstellar wind model to describe the radio 
light curve for SN 1998bw. This is in many ways similar to the Chevalier jl4j 
model for Type Ib/c SNe, that also included synchrotron self-absorption (SSA). 
To produce the first bump in the radio light curves of SN1998bw, Li & Chevalier 
postulated a boost of blastwave energy by a factor of ~ 2.8 on day ~22 in the 
observer's time frame. They did not discuss the second bump. 

Modeling of the radio data for S N 199 8bw with the well established formalism 
for RSNe presented by Weiler et al.|126] shows that such an energy boost is not 
needed. A fast blastwave interacting with a dense, slow, stellar wind-established, 
ionized CSM, that is modulated in density over time scales similar to those seen 
for RSNe, can produce a superior fit to the data. No blastwave re-acceleration is 
required and no SSA at early times is apparent. The parameters of the best fit 
model are given in Table [l] and shown as the curves in Fig. [l| A visual comparison 
of the curves in Fig. |l| with those of Li & Chevalier |35| Fig. 9, shows that the 
purely thermal absorption model with structured CSM provides a superior fit. 



Table 1. SN1998bw/GRB980425 Modeling Results 1 

Parameter Value 

q (spectral index) —0.71 

P (decline rate) —1.38 

Kl 2.4 x 10 3 

K 2 
5 

Kl 1.7 x 10 3 

5' -2.80 

K 4 1.2 x 10" 2 
to (Explosion Date) 1998 Apr. 25.90915 

(teem peak - to)(days) 13.3 

S6cm poak(mjy) 37.4 

d(Mpc) 38.9 

L 6 cm peak (ergs s" 1 Hz" 1 ) 6.7 x 10 28 

M(M Q yr" 1 ) 6 2.6 x 10 _s 



1 See Weiler et al. [ 126 for parameter definitions. 

"Enhanced by a factors of 1.6 {K^) and 2.0 (-K3), corresponding to a density increase 
of 40%, over the intervals day 25 - 75 and day 165 - 215, although the latter interval 
could be as long as 100 days and still be compatible with the available data. 
'Assuming U = 23 days, t — (teem peak — to) = 13.3 days, m — — (a — f3 — 3)/3 = 0.78, 
«Jwind = 10 km s" 1 , Vi = ^blastwave = 230, 000 km s _1 , T = 20, 000 K, average number 
of clu mps along the line-of-sight iV = 0.5 and volume filling factor cj> — 0.22 (see Weiler 
et al.JLM], Case 3). 
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One should note that the fit listed in Table [l] and shown as the curves in Fig. [I] 
requires no "uniform" absorption (K2 — 0) so all of the free-free (f-f ) absorption 
is due to a clumpy medium as described in Equations 1 and 5 of [126|. These 
results, combined with the estimate of a high blastwave velocity, suggest that the 
CSM around SN1998bw is highly structured with little, if any, inter-clump gas. 
The clump filling factor must be high enough to intercept a considerable fraction 
of the blastwave energy and low enough to let radiation escape from any given 
clump without being appr ecia bly absorbed by any other clump which is Case 3 
discussed by Weiler et al.|126|. The blastwave can then easily move at a speed 
that is a significant fraction of the speed of light, because it is moving in a very 
low density medium, but still cause strong energy dissipation and relativistic 
electron acceler ation at the clump surfaces facing the SN explosion center. 

Weiler et al .j 126 also noted from the fit given in Table |l] that the presence of 
a K4 (see |l26| ) factor implies there is a more distant, uniform screen of ionized 
gas surrounding the exploding system that is too far to be affected by the rapidly 
expanding blastwave and provides a time independent absorption. 



1.5 Physical Parameter Estimates 



Using the fitting parameters from Table [l] and Equations 11 and 16 of [126] 



Weiler et al. [126 1 estimated a mass-loss rate from the preexplosion star. The 
proper parameter assumptions are rather uncertain for these enigmatic objects 
but, for a preliminary estimate, they assumed t\ — 23 days, t — (te cm p0 ak — to) = 
13.3 days, m = -(a - (3 - 3)/3 = 0.78 (Equation 7 of p2^), w wind = 10 
km s _1 (for an assumed RSG progenitor), V\ — Wbiastwave = 230,000 km s _1 , 
and T = 20,000 K. They also assumed, because the radio emission implies that 
the CSM is highly clumped (ie K2 — 0), that the CSM volume is only sparsely 



occupied (N = 0.5, <j> = 0.22; see Weiler et al.|26|, Case 3). Within these 
rather uncertain assumptions, Equations 11 and 16 of Weiler et al.[ [126| yield an 
estimated mass-loss rate of M ~ 2.6 x 10~ 5 M Q yr -1 with density enhancements 
of ~ 40% during the two known, extended "bump" periods. 

Assuming that the blastwave is traveling at a constant speed of ~ 230,000 
km s _1 the beginning of the first "bump" on day 25 implies that it starts at 
~ 5.0 x 10 16 cm and ends on day 75 at ~ 1.5 x 10 17 cm from the star. Cor- 
respondingly, if it was established by a 10 km s _1 RSG wind, the 50 days of 
enhanced mass- loss ended ~ 1, 600 yr and started ~ 4, 700 yr before the explo- 
sion. The earlier high mass-loss rate epoch indicated by the enhanced emission 
on day 192 in the measurement gap between day 149 and day 249 implies, with 
the same assumptions, that it occurred in the interval between ~ 9, 400 yr and 
~ 15, 700 yr before explosion. It is interesting to note that the time between the 
presumed centers of the first and second increased mass-loss episodes of ~ 9, 400 
yr is comparable to the ~ 12, 000 yr before explosion at which SN1979C had a 
significant mass-loss rate increase |7(J and SN1980K had a significant mass-loss 
rate decrease |7j|, thus establishing a possible characteristic time scale of ~ 10 4 
yr for significant changes in mass-loss rate for preexplosion massive stars. 
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1.6 Radio Emission from SN1998bw/GRB980425 and Other GRB 
Radio Afterglows 

Because the suggestion of a possible relation between SN1998bw and GRB 
980425, it has remained a tantalizing possibility that the origin of at least some 
GRBs is in the better known Type Ib/c SN phenomenon. First, of course, one 
must keep in mind that there may be (and probably is) more than one origin 
for GRBs, a situation that is true for most other classes of objects. For example, 
SNe, after having been identified as a new phenomenon in the early part of the 
last century, were quickly split into several subgroups such as Zwicky's Types I 
- V, then coalesced back into just two subgroupings based on Ha absent (Type 
I) or Ha present (Type II) in their optical spectra. This simplification has not 
withstood the test of time, however, and subgroupings of Type la, lb, Ic, II, lib, 
Ilpec, Iln, and others have come into use over the past 20 years. 

GRBs, although at a much earlier stage of understanding, have similarly 
started to split into subgroupings. The two currently accepted groupings are 
referred to as "fast-hard" and "slow-soft" from the tendency of the 7-ray emission 
for some to evolve more rapidly (mean duration ~0.2 s) and to have a somewhat 
harder spectrum than for others that evolve more slowly (mean duration ~20 s) 
with a somewhat softer spectrum p4fl . 

Because we are only concerned with the radio afterglows of GRBs here, all 
of our examples fall into the slow-soft classification, at least partly because the 
fast-hard GRBs fade too quickly for followup observations to obtain the pre- 
cise positional information needed for identification at longer wavelengths. It 
is therefore uncertain whether fast-hard GRBs have radio afterglows or even 
whether the slow-soft GRBs represent a single phenomenon. If, however, we as- 
sume that at least some types of slow-soft GRBs have a similar origin and that 
GRB980425/SN1998bw is a key to this puzzle telling us that such "slow-soft" 
GRBs have their origin in at least some types of SNe, we can investigate relations 
between the two observational phenomena. 



2 Gamma-Ray Bursters 
2.1 Radio Detections 

Beyond GRB980425, there are relatively few GRBs with detected radio after- 
glows and only six of these as of 2001 Dec. 31 have sufficient radio light curve 
information to permit approximate model fits. Additionally, because the opti- 
cal/infrared (OIR) data appear consistent with a synchrotron origin similar to 
that of the radio emission, we have collected the available OIR data and included 
it in our model fitting to better constrain the source parameters. Although de- 
tailed modeling is beyond the scope of this review, we apply the parameterization 



of Equations 1 to 10 of Weiler et al.[126| to the available radio and OIR data 



in attempt to highlight some of the gross properties of the GRB afterglow pro- 
cesses. Because the OIR data suffer an extinction that is absent in the radio, 
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a zero redshift color excess [E(B — V)] was also obtained from the fitting by 
adopting the Galactic law of Savage & Mathis |97j . 

GRB970508 was discovered by the Beppo-Sax team on 1997 May 8.904 
UT . These results showed detection of an afterglow in all wavelength bands 
including X-ray p6| |, optical Q, and radio p5fl . Frail & Kulkarni |2f| derived a 
position from their 8.46 GHz VLA observations of RA(J2000) = 06*53" l 49!45, 
Dec(J2000) = +79°16T9'/5 with an error of -O'.'l in both coordinates. Metzger 
et al.[Q found a redshift of z = 0.835, which Bloom et al.{| confirmed for the 
host galaxy. 

The radio data were obtained from the references [3 30,37.101 ,[U)5i and the 
OIR data were obtained from the references §,|l|,|l|j37jt ' 



33,[41|,£6|]9J 



108 j . Repre- 



sentative data for GRB970508 are plotted, along with curves from the best fit 
model, in Fig. || and the parameters of the fit are listed in Table ||. 



Table 2. Fitting Parameters for GRB Radio Afterglows 1 
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GRB970508 


0.08 


0.00 


0.835 


-0.63 


-1.18 


1.35 x 


w' 1 


1.81 x 


10 3 


-1.75 


1.39 


x 10 31 


GRB980329 


0.07 


0.28 


= 1.000 


-1.33 


-1.09 


1.54 x 


10 4 


1.33 x 


10 5 


-1.16 


7.69 


x 10 30 


GRB980425 a 






0.0085 


-0.71 


-1.38 


2.37 x 


10 3 


1.73 x 


10 3 


-2.80 


6.70 


x 10 2S 


GRB980519 


0.27 


0.00 


= 1.000 


-0.75 


-2.08 


8.45 x 


10 1 


1.37 x 


10 1 


-3.57 


7.52 


x 10 30 


GRB991208 


0.02 


0.07 


0.707 


-0.58 


-2.27 


5.11 x 


10 2 


3.53 x 


10 3 


-3.29 


2.07 


x 10 31 


GRB991216 


0.63 


0.11 


1.020 


-0.28 


-1.38 


7.07 x 


10° 


2.50 x 


K) 1 


-1.40 


1.05 


x 10 31 


GRB000301C'' 


0.05 


0.05 


2.034 


-0.60 


= -1.75 


2.31 x 


10 2 


2.72 x 


10 3 


-2.06 


3.77 


x 10 31 



1 The fits do not generally require use of K2, 8, K4, K$, 8" , Ke, 8"' (see Weiler et al. 126]) 
so that these parameters are not included. 

2 Gala ctic extinction in the direction of the 7-ray burster (GRB) was obtained from 

3 A . 



Additional extinction, at the host galaxy redshift (see §2.1), is needed in some cases 
to provide a good fit to the optical/infrared (OIR) data. 
4 Where unknown, the redshift is defined to be z = 1.000. 
5 Derived for 6 cm in the rest frame of the observer, not the emitter. 
"The best fit includes a K 4 = 1.24 x 10" 2 term. 

6 The OIR data and radio data appear to have different rates of decline (/3oir ~ —2.0, 
/Sradio ~ —1.5) implying that there may be a break between the two regimes. For the 
purposes of this review the average of /3 = —1.75 has been adopted. 



Examination of Figs. §a,b shows that the parameterization listed in Table || 
describes the data well in spite of the very large frequency and time range. In 
Fig. ||a the 232 GHz upper limits and the 86.7 GHz limits and detections are in 
rough agreement with the model fitting; the 8.5 GHz and 4.9 GHz measurements 
are described very well, and even the 1.5 GHz data are consistent with the 
parameterization if significant interstellar scintillation (ISS) is present (see § 
[2.3| ). Waxman et al. 120 have already ascribed the large fluctuations in the flux 
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Fig. 2. GRB970508 at radio wavelengths (left, 2a) of 1.3 mm (232 GHz; open circles, 
solid line), 3.5 mm (86.7 GHz; crosses, dashed line), 3.5 cm (8.5 GHz; open squares, 
dash-dot line), 6 cm (4.9 GHz; open triangles, dotted line), and 20 cm (1.5 GHz; open 
diamonds, dash-triple dot line) and optical/IR bands (right, 2b) of B (681 THz; circles, 
solid line), V (545 THz; crosses, dashed line), R (428 THz; open squares, dash-dot line), 
I (333 THz; open triangles, dotted line), and K (138 THz; open diamonds, dash-triple 
dot line). (NB: To enhance clarity, the scales are different between the two figures even 
though the fitting parameters are the same and not all available bands are plotted even 
though they were used in the fitting.) 



density at both 8.5 and 4.9 GHz to ISS, and one expects such ISS to also be 
present at 1.4 GHz. 

In Fig. ||b, although the OIR data show significant scatter at individual 
frequencies, the data are consistent with a non-thermal, synchrotron origin that 
has the same decline rate (/?) and spectral index (a) as in the radio regime, 
indicating that no spectral breaks have occurred between the two observing 
ranges. The color excess of E(B — V) — 0.08 mag from Schlegel et al.[10C] is 
consistent with the best fit to the data, implying little extinction in the host 
galaxy. 

The most obvious characteristic of the OIR data is that in the initial interval 
between the time of the GRB and day 1.75 the data are not well modeled by 
the parameterization. This is not surprising, because the modeling contains no 
parameters to describe any turn-on effects for the synchrotron emission. Thus, 
the initial OIR data prior to day 1.75 are particularly suitable for constraining 



explosion models such as those discussed in § 2.4. 

GRB980329 was discovered by the Beppo-Sax team (BATSE Trigger # 
6665, @) on 1998 Mar. 29.1559 UTjp). The afterglow was detected in all 
wavelength bands including X-ray [ 13 1 1 , optical [ pi] , and radio ]110| . Taylor 
et al.[ |110| derived a position from their VLA observations of RA(J2000) = 
07' l 02 m 38?022, Dec(J2000) = +38°50'44'.'02 with an uncertainty of ±0'.'05 in 
each coordinate. Unfortunately, no redshift has been obtained for the optical 
afterglow of GRB980329 or its inferred parent galaxy so that, except for argu- 
ments that is quite distant with, perhaps, z ~ 5, no reliable distance estimate is 
available. 
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Fig. 3. GRB980329 at radio wavelengths (left, 3a) of 0.9 mm (352 GHz; open circles, 
solid line), 3.5 cm (8.5 GHz; crosses, dashed line), 6 cm (4.9 GHz; open squares, dash- 
dot line), and 20 cm (1.4 GHz; open triangles, dotted line) and optical/IR bands (right, 
3b) of R (428 THz; open circles, solid line), I (333 THz; crosses, dashed line), J (240 
THz; open squares, dash-dot line), and K (136 THz; open triangles, dotted line). (NB: 
To enhance clarity, the scales are different between the two figures even though the 
fitting parameters are the same and not all available bands are plotted even though 
they were used in the fitting.) 



The radio data were obtained from [ |104|]l0ii ,111] and the OIR data were 
obtained from pB| , po| . Representative data for GRB980329 are plotted, along 
with curves from the best fit model, in Fig. [| and the parameters of the fit are 
listed in Table §. 

Examination of Figs. |^a,b shows that the parameterization listed in Table || 
describes the data rather well over the broad parameter space in time and fre- 
quency. In Fig. ||a the 352 GHz detections and upper limits are in good agree- 
ment with the parameterization; the 8.5 GHz and 4.9 GHz measurements are 



described very well although there may be some ISS present (see § 2J3); the 1.4 
GHz upper limits are consistent with the parameterization. 

In Fig. ||b the data are mostly upper limits, with only a few detections at 
K-band (136 THz) and R-band (428 THz). However, these are surprisingly well 
described by the parameterization if a zero redshift color excess of E(B — V) = 
0.75 mag is assumed. This value for color excess is much higher than the value 
of E(B — V) = 0.073 mag from Reichart et al.|M| for Galactic color excess in 
that direction. However, if we assume that the additional color excess arises in 
the host galaxy and that the color excess scales with redshift as (1 + z)~ 1,25 , as 
is appropriate for an adopted Small Magellanic Cloud extinction law |87j, then 
the host galaxy contributes E(B — V) = 0.28 mag (Table ||). Such a color excess 
is fairly typical for t he disk of a late-type galaxy. 

Taylor et al.[lll have also modeled the radio data and invoked a somewhat 
steeper, inverted spectrum with a = — 1.7 between 4.9 and 8.3 GHz, flattening to 
a = —0.8 between 15 and 90 GHz caused by a SSA component with a turnover 
frequency near 13 GHz. Extrapolating to higher frequencies, their model predicts 
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Fig. 4. GRB980519 at radio wavelengths (left, 4a) of 3.5 cm (8.5 GHz; open circles, 
solid line), 6 cm (4.9 GHz; crosses, dashed line), and 20 cm (1.4 GHz; open squares, 
dash-dot line) and optical/IR bands (right, 4b) of B (681 THz; open circles, solid line), 
V (545 THz; crosses, dashed line), R (428 THz; open squares, dash-dot line), and I 
(333 THz; open triangles, dotted line). (NB: To enhance clarity, the scales are different 
between the two figures even though the fitting parameters are the same and not all 
available bands are plotted even though they were used in the fitting.) 



a rather low 350 GHz flux density of only ~1.7 mJy that is incompatible with 
the James Clerk Maxwell Telescope measurements. 

Such complexity is not be needed, however. Our parameterization listed in 
Table |^ and shown in Fig. |^a yields a good description of the data and pre- 
dicts a 350 GHz flux density of ^3.0 mJy, in much better agreement with the 
observations. 

Our model fit to the available first 30 days of radio data indicates that 
GRB980329 should be detectable with the VLA at centimeter wavelengths, with 
little decline, for an extended period. Although prediction of exact flux densities 
at later times is not reliable because the decline phase of the model is not well 
constrained by the few available data, it is interesting to note that GRB980329 
has apparently been detected by the V LA at 8.46 GHz, 4.86 GHz, and 1.43 GHz 
for up to 500 days after the outburst [ 130 1 . 

GRB980519 was discovered by the Beppo-Sax team (BATSE Trigger # 
6764) on 1998 May 19.51410 UT @. The afterglow was detected in all wave- 
length bands including X-ray J78|, optical Q, and radio @. Frail et al.|p6| 
derived a position from their VLA observations of RA(J2000) = 23' l 22 m 2i?49, 
Dec(J2000) = +77°15'43'.'2 with an uncertainty of ±0'.'l in each coordinate. Un- 
fortunately, no redshift has been obtained for the optical afterglow so that no 
distance estimate is available. 



The radio data were obtained from |3C[ and the OIR data from |46] , |55 ,117| 
Representative data for GRB980519 are plotted, along with curves from the best 
fit model, in Fig. |] and the parameters of the fit are listed in Table ||. 

Examination of Figs. ||a,b shows that the parameterization listed in Table || 
describes the data reasonably well, even though the radio data shown in Fig. 
0a are very limited and of relatively poor quality. For example, the reported 
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Fig. 5. GRB991208 at radio wavelengths (left, 5a) of 1.2 mm (250 GHz; open circles, 
solid line), 2.0 cm (15.0 GHz; crosses, dashed line), 3.5 cm (8.5 GHz; open squares, 
dash-dot line), 6 cm (4.9 GHz; open triangles, dotted line), and 20 cm (1.4 GHz; open 
diamonds, dash-triple dot line) and optical/IR bands (right, 5b) of B (681 THz; open 
circles, solid line), V (545 THz; crosses, dashed line), R (428 THz; open squares, dash- 
dot line), and I (333 THz; open triangles, dotted line). (NB: To enhance clarity, the 
scales are different between the two figures even though the fitting parameters are the 
same and not all available bands are plotted even though they were used in the fitting.) 



detections at 4.9 GHz after day 50 are barely more than 3<r and thus of limited 
reliability, so that the fit is not well constrained and the significance of their 
deviation from the best fit curve is unknown. The data at both 8.5 GHz and 
4.9 GHz have significant fluctuations yielding detections and 3tr upper limits at 
relatively small time separations so that GRB980519 may be undergoing ISS 
(see § 2.3). Only upper limits are available at 1.4 GHz, but they are consistent 
with the best fit model. 

In Fig. |Jb there is good coverage in both frequency and time. Although there 
is some fluctuation, the data are rather well described by the parameterization 
if a color excess of E(B — V)= 0.25 mag is assumed. This is r ather close to the 
value of E(B — V) — 0.267 mag suggested by Schlegel et al. |l00| in that direction 
from their IRAS 100 ^m maps. 

GRB991208 was discovered by the Ulysses and KONUS and NEAR teams 
on 1999 Dec. 08.1923 UT @. The afterglow was detected in the optical § 
and radio [^7] wavelength bands. There does not appear to have been an X- 
ray detection. Frail & Kulkarni p7| derived an position for the radio transient 
of RA(J2000) = 16 h 33 m 53'/50, Dec(J2000) = +46°27'20"9 with no error given, 
but normally for the VLA positional errors are < O'.'l. Dodonov et al.|22| found 
a redshift for the parent galaxy of z = 0.707 ± 0.002. 

The radio data were obtained from Galama et al.|Q and the OIR data from 
Castro- Tirado et al.jljj]. Representative data for GRB991208 are plotted, along 
with curves from the best fit model, in Fig. ^| and the parameters of the fit are 
listed in Table §. 

Examination of Figs. |^a,b shows that the parameterization listed in Table ^| 
describes the data reasonably well. In Fig. ||a the radio data are quite limited 
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and have relatively large scatter. For example, the reported detection at 15 GHz 
on day 7.4 is preceded by a 3<7 upper limit of comparable magnitude on day 5.4 
and followed by an upper limit on day 12.5, so that the reality of the day 7.4 
detection must be called into question. In any case, although the fit is not well 
constrained by the radio data, the model curves describe its evolution reasonably 
well. The data are too sparse to judge if ISS is present. 

In Fig. ||b there is good coverage in both frequency and time. There is some 
fluctuation, particularly for the four R-band observations on day 4. However, 
examination of the data reveals that the three high R-band measurements on 
day 4.026, 4.058, and 4.096 have significantly larger errors (<~ 0.4 mag) than the 
R-band datum on day 4.068 (error = 0.050 mag) that is consistent with the fitted 
curves. In general, the data are rather well described by the parameterization 
if a zero redshift color excess of E(B — V) = 0.15 mag is assumed. This is 
appre ciab ly higher than the value of E(B— V) = 0.016 mag suggested by Schlcgcl 
et al.[10C] in that direction from their IRAS 100 /i maps, implying, under the 
same assumptions as were used for GRB980329 above, a E(B — V) = 0.07 mag 
at the redshift of the host galaxy (Table ||) . 

GRB991216 was discovered by BATSE on 1999 Dec. 16.671544 UT (BATSE 
trigger # 7906) The afterglow was detected in all wavelength bands includ- 
ing x-ray |109| , optical [115|, and radio [ 1 1 2 1 . Taylor & Berger |112| derived a 
position from their 8.5 GHz VLA observations of RA(J2000) = 05 7I 09 m 31!297, 
Dec(J2000) = +11°17'07'.'25 with an error of ±0'.'l in each coordinate. Vreeswijk 
et ai |118| suggested a redshift of z > 1.02 based on the highest redshift of three 
possible absorption systems seen in the optical afterglow. 

The radio data were obtained from Frail et al.lBll and the OIR data from 



Halpern et al.[]47|, Frail et al. 31 , and Garnavich et al.|4j. Representative data 
for GRB991216 are plotted, along with curves from the best fit model, in Fig. || 
and the parameters of the fit are listed in Table |^. 

Examination of Figures |^a,b shows that the parameterization listed in Table ^| 
fits the data reasonably well, even with the radio data (Fig. |^a) being quite 
limited and with relatively large scatter. The single detection at 8.5 GHz on 
day 42.49 surrounded by much lower 3cr upper limits on days 38.28 and 50.51, if 
correct, is not well described. The data are too sparse to judge if ISS is observed. 

In Fig. ^|b there is reasonably good coverage in both frequency and time, 
particularly at R-band. A satisfactory fit to the data requires a rather high zero 
redshift color excess of E(B — V) = 0.90 mag. This is significantly hig her than 
the value of E(B — V) = 0.63 mag suggested by Schlegel et al. |l00| in that 
direction from their IRAS 100 fi maps or the value of E(B — V) = 0.40 mag 
derived by Halpern et al.[[l7| from the Galactic 21 cm column density in that 
direction. If, as above, we assume that the additional extinction arises in the 
host galaxy, then with the same assumptions as were used for GRB980329, we 
obtain a color excess of E(B — V) =0.11 mag at the redshift of the host galaxy 
(Table |). 

GRB000301C was discovered by the Ulysses and NEAR teams on 2000 
Mar. 01.4108 UT pi. The afterglow was detected in the optical p| and radio 
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(t-y/days [t B = (IB-Dec-99) +0.000 days.] (t-t )/day5 [t a = (16-Dec-99) +0.000 days.] 



Fig. 6. GRB991216 at radio wavelengths (left, 6a) of 0.9 mm (350 GHz; open circles, 
solid line), 2.0 cm (15.0 GHz; crosses, dashed line), 3.5 cm (8.5 GHz; open squares, 
dash-dot line), 6 cm (4.9 GHz; open triangles, dotted line), and 20 cm (1.4 GHz; open 
diamonds, dash-triple dot line), and optical/IR bands (right, 6b) of R (428 THz; open 
circles, solid line), I (333 THz; crosses, dashed line), J (239 THz; open squares, dash-dot 
line), and K (136 THz; open triangles, dotted line). (NB: To enhance clarity, the scales 
are different between the two figures even though the fitting parameters are the same 
and not all available bands are plotted even though they were used in the fitting.) 



H wavelength bands. There does not appear to have been an X-ray detection. 
Fynbo et al.(33| derived a position for the optical transient of RA(J2000) = 
16' l 20 m 18'.'6, Dcc(J2000) = +29°26'36" with an error of - 1" in both coordi- 
nates. Castro et al.jl(| (see also [ L02 1) found a redshift for the parent galaxy of 
z = 2.0335 ±0.003. 

The radio data were obtained from Berger et al.Q and the OIR data from 
Sagar et al.[p5), Bhargavi & Cowsik ||, Masetti et al.|39|, Rhoads & Fruchter 
[ j92| , and Jensen et al.pq]. Representative data for GRB000301C are plotted, 
along with curves from the best fit model, in Fig. |?] and the parameters of the 
fit are listed in Table ||. 

Examination of Figures 0a,b shows that the parameterization listed in Ta- 
ble H is successful in that it describes the data rather well over the large frequency 
and time range. In Fig. |a the 350 GHz upper limits and the 250 GHz detec- 
tions and limits are in rough agreement with the parameterization; the 15 GHz 
and 8.5 GHz measurements are also consistent with the modeling, although the 
model drops off a bit faster at late times than the data and the fit would be 
improved by a somewhat flatter decline rate of (3 ~ —1.5. The 1.4 GHz data are 
generally consistent with the model, although ISS (see § 2.3) may be present in 
the first measurement. 

In Fig. ^b, although the OIR data show significant scatter at individual 
frequencies, they are consistent with a non-thermal, synchrotron origin. A decline 
rate of (3 ~ —2.0, somewhat steeper than in the radio, would provide an improved 
fit. The fit indicates a best value for the zero redshift color excess of E(B — V) = 
0.25 mag, significantly higher than the value of E(B — V) — 0.052 mag suggested 
by Schlegel et al. 1 100 in that direction from their IRAS 100 /i maps. As above, 




Fig. 7. 000301C at radio wavelengths (left, 7a) of 0.9 mm (350 GHz; open circles, 
solid line), 1.2 mm (250 GHz; crosses, dashed line), 2.0 cm (15.0 GHz; open squares, 
dash-dot line), 3.5 cm (8.5 GHz; open triangles, dotted line), and 6.0 cm (4.9 GHz; open 
diamonds, dash-triple dot line) and optical/IR bands (right, 7b) of U (832 THz; circles, 
solid line), B (681 THz; crosses, dashed line), R (428 THz; open squares, dash-dot line), 
I (333 THz; open triangles, dotted line), and K (136 THz; open diamonds, dash-triple 
dot line). (NB: To enhance clarity, the scales are different between the two figures even 
though the fitting parameters are the same and not all available bands are plotted even 
though they were used in the fitting.) 



this implies a color excess of E(B — V) — 0.05 mag at the redshift of the host 
galaxy (Table |). 

The most obvious characteristic of the OIR data is that the initial interval 
between the time of the GRB and day 3.5, is not well modeled by the parame- 
terization. This is not surprising, because the modeling contains no parameters 
to describe any turn-on effects for the synchrotron emission. A similar turn-on 
effect was seen for the early OIR data for GRB970508, although not for the other 
GRBs discussed here. 

It should be noted that the possible steepening of the flux density decline 
rate from j3 ~ —1.5 (radio) to f3 ~ —2.0 (OIR) may indicate a break in the 
decline rate somewhere between the two observing ranges, possibly owing to 
synchrotron losses of the high energy electrons. 

2.2 Parameterization 

Even though the study of the GRB phenomenon is still very much in its early 
stages and, as has been pointed out above, the GRBs for which we have any 
afterglow information are all of the "slow-soft" category with those about which 
we have significant radio information only a small subset, the parameterization 
studies (see Table ||) allow us to draw a some tentative conclusions: 

1 . Of the seven relatively well observed radio GRBs (RGRBs) (including GRB 
980425/SN1998bw), five appear to have similar relativistic electron acceler- 
ation processes that generate a spectral index a ~ —0.6 to ~ —0.7. The two 
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possible exceptions are GRB980329 (a 1.3) and GRB991216 (a 0.3) 

which, interestingly enough, straddle the average spectral index of the other 
objects. The average spectral index value is very similar to the canonical 
value for extragalactic radio sources. 

The RGRBs can be described rather well from radio through OIR by the 
same formalism of Equations 1 to 10 of Weiler et al. [126 that is successful in 
describing the radio emission from RSNe as arising from a blastwave/CSM 
interaction. The obvious conclusion is that at least this subset of GRBs arise 
from the explosion of a compact, presumably stellar-sized object embedded 
in a dense CSM. 

Only two of the objects, GRB970508 and GRB000301C, appear to have been 
observed in a phase for which they show the turn-on of their OIR radiation. 
Such a turn-on is presumably related to the unmodeled relativistic electron 
acceleration at very early times and provides a test of explosion models. 
None of the objects requires the inclusion of a uniformly distributed absorb- 
ing component to the early time radio absorption (K 2 in Equation 1 of |l26| ) 
that implies a highly filamentary or clumpy structure in the CSM with which 
the blastwave is interacting. 

Estimates for the color excess, E(B — V), obtained from the modeling imply 
that many of the objects suffer appreciable extinction at the redshift of 
their parent galaxies because fitted zero redshift values significantly exceed 
estimates of Galactic extinction in their directions. The values for the excess 
extinction at the redshifts of the host galaxies noted above and in Table || 
are rather typical for late-type galaxies. 

All of the RGRBs except for GRB980425/SN1998bw have similar observed 
peak 6 cm radio spectral luminosities of ~ 10 31 ergs s _1 Hz -1 . This implies 
that GRBs are rare objects for which one must survey very large volumes of 
space in order to detect examples. Such examples are, therefore, generally 
at very large distances and we detect only the brightest members of the 
luminosity distribution. GRB980425/SN1998bw is the only known nearby, 
lower luminosity exception. 

The prevalence of interstellar scintillation (ISS) in the early radio emission 
from many of these objects implies that the RGRBs are initially very small 
in angular size although they expand rapidly. 



2.3 Discussion 

Although many useful conclusions can be drawn from the light curve parame- 
terizations, one should keep in mind that a more detailed study of radio GRBs 
will have to take into account several physical effects that are not seen in RSNe, 
and that have not been included in this brief overview, but need further consid- 
eration. 



Interstellar Scintillation: Because their high radio luminosity and low ab- 
sorption allows detection at great distance and quite early when they are still 
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of very small angular size, GRB radio afterglows appear to be so compact as to 
exhibit ISS during the first few days or weeks of detectability. This possibility 
was proposed by Goodman Q for GRB970508 based on earlier work by Rickett 
[ |9l| for pulsars. After consideration of the several regimes of strong, weak, re- 
fractive, and diffractive scattering, Goodman [Q concluded that both diffractive 
and refractive scintillation are possible for radio afterglows and that observation 
of the effects of scintillation can place limits on their /Ltas (micro-arcsecond) sizes 
at levels far too small to be resolved with VLBI. 

GRB970508 shows strong flux density fluctuations at bot h 8.4 6 and 4.86 



GHz until age ~4 weeks after explosion, that Waxman et al. [ 120 1 attributed 
to diffractive scintillation. After ^1 month, they found that the modulation 
amplitude decreased, which is consistent with the diffractive scintillation being 
quenched by the increased size of the radio emitting region. They also took this 
increasing source size to be consistent with, and supportive of, the "fireball" 
mode l pre dictions (see § From their 4.86 and 8.46 GHz results, Waxman 



et al. [ 12C | concluded that the quenching of diffractive scintillation at ^4 weeks 
implies a size at that epoch of ~ 10 cm and an expansion speed comparable 
to that of light. 



In contrast to the conclusion of Waxman et al. [ 120 that GRB970508 is un- 



dergoing strong diffractive scintillation during the first 30 days after explosion, 



Smirnova & Shishov 1 103 ] concluded that the radio afterglow is, in fact, under- 
going only weak scintillation at 4.86 and 8.46 GHz at early times with refractive 
scintillation dominating at 1.43 GHz. 

GRB980329 shows rapid flux density fluctuations at 4.9 and 8.3 GHz that 
are extinguished by age ~3 weeks. Although they did not analyze the scintillation 



data in detail, Taylor et al.|lll] pointed out a similarity to the better studied 
scintillations of GRB970508 and suggested that the early-time angular size of 
GRB980329 may be even smaller than the ~ 3 /Ltas inferred for GRB970508 by 



Goodman E4] and Waxman et al.[12C| because GRB980329 is at a lower Galactic 
latitude, that should more quickly quench ISS. 

GRB980519 shows strong modulation of the flux density at 4.86 and 8.46 
GHz during the first ~20 days after the GRB burst, that Frail et al.[^o) in- 
terpreted as being due to diffractive ISS. They derived a resulting maximum 
radio source size of < 0.4 /ias, an extremely compact object. As was seen for 
GRB970508, the 1.4 GHz emission from GRB980519 seems to be suppressed, 
in this case below their detectability limit because their three measurements at 
1.4 GHz are all upper limits. Frail et al.|^(| attributed this suppression to syn- 
chrotron selfabsorption (SSA) with a turn-over frequency between 1.43 and 4.86 
GHz. 



Cosmological Effects: Because the GRBs are at cosmological distances (z ~ 
1), there are two effects that must be taken into account: 

1. There is a time dilation that slows the light curve evolution in the observer's 
frame with respect to the time evolution in the emission frame. This is a 
straightforward correction and has been elaborated by a number of authors 
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(see, e.g., p"7| , [48| , |33| ) . The time dilation results in a true emitted time to 6 
cm peak flux density [(i6cm peak — £o)cmit] shorter than that actually observed 
[(i6cm peak — io)obs]- The correction takes the form 

(^6cm peak ^o)cmit — (^6cm peak ^o)obs (1 ~\~ (1) 



and must be applied to the measured times from explosion to 6 cm peak flux 
density to obtain true times. 
2. There is a correction of the observed flux density, (S l 6 C m)obs, owing to the red- 
shift making the observed frequency different from that emitted, (S , 6 C m)emit) 
for sources with non-zero spectral indices. The correction normally takes the 
form, for S oc v +a , 

(5*60111)611111 — (S 6cm Us (1 + z)- a (2) 

However, Chevalier & Li jl6) have proposed an "equality of peaks" on the- 
oretical grounds so that such a correction may be less important than ex- 
pected. 



Relativistic Effects: Essentially all researchers agree that the GRB phe- 
nomenon involves relativistic motion of the emitting region. Whether the motion 
involves a spherical, relativistic fireball or a directed relativistic jet is probably 
not of concern from the radio observer's standpoint because the emission from 
the CSM interaction is probably not highly directed. 

For relativistic corrections there are two factors - the relativistic motion 
factor r [r = (1 — and, if the motion is directed, the viewing angle 9. If, 

as expected, the radio emission is not highly directed, then taking 8 — leaves 
only the relativistic motion factor r to be determined. 

The factor r is difficult to estimate from observations. Theoretical modeling 
predicts F of several hundred very early in the expansion phase Jl 19[ declining 
to subrelativistic motion after only a few days to a few weeks p9|j50| . However, 
although highly speculative, there is, perhaps, one method for estimating r. 

Weiler et al.|126] pointed out that by assuming a relation between Type lb 
SN1998bw and GRB980425, they could obtain estimates of T - 1.6 - 2.0 for 
the radio afterglow. They obtained one estimate by postulating that SN1998bw 
was expanding at ^230,000 km s if its brightness temperature at 6 cm peak 
flux density was 10 K, as required by the inverse Compton limit. The second 
estimate was obtained by assuming that the true 6 cm peak spectral luminosity 
was the same as that for the average of all known Type Ib/c RSNc. Then, 
if one assumes: (a) all GRB radi o afterglows arise in Type Ib/c SNe and (b) 
the postulate of Weiler et al. |l25| that all Type Ib/c SNe may be approximate 
standard radio candles with an Lq ciu peak ~ 1-3 x 10 27 erg s _1 Hz -1 , and (c) any 
increase in the observed flux density is solely due to relativistic boosting, then 
one can estimate a relativistic r factor for each GRB radio afterglow. 

Applying these assumptions to the peak 6 cm radio luminosity (L ecm pea k) 
values obtained for the GRBs and listed in Table ||, we obtain r values ~ 10. 
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Such fs arc similar to those found for VLBI "super luminal" sources (see e.g., 



9p9fl) and for the Galactic micro-QSO superluminal sources p^,p4|,[74[ll3] 



2.4 Theoretical Models 

The origin of 7-ray bursts is a long standing problem dating back to the 1960s. 
With little actual data to constrain them, many theories were developed. Rud- 
erman [^3) first summarized these in the mid-1970s. Cavallo & Rees were 
probably the first to use the term "fireball" in an article discussing several pos- 
sibilities for the conversion of massive amounts of injection energy into 7-rays, 
including the possibility of relativistic expansion and the conversion of kinetic 
energy into 7-ray luminosity through the impact on an external medium. 

The 1980s brought the more detailed calculations of Paczynski J79| and Good- 
man Q that suggested that the GRBs might be at cosmological distances of 
z ~ 1 or more and that the energies involved were at least comparable to that 
expected from an SN (~ 10 51 ergs). In the 1990s, of course, the field finally came 
of age with the Compton Gamma-Ray Observatory to provide almost daily de- 
tections of GRBs and Beppo-Sax to provide rapid followup for many bursts and 
accurate position determination for X-ray afterglows. This increased data flow 
refined the model studies for the origin of GRBs so that there are basically 
only two models presently under serious consideration - relativistic fireball and 
relativistic jet. Because of the observational (particularly radio) orientation of 
this review, we shall not discuss the mechanisms for the source of the energy 
release, whether neutron star or black hole formation in SNe, coalescence of 
compact objects, or more exotic phenomena. We shall be more concerned with 
the observable consequences of this great energy release and, where possible, 
similarities to known phenomena. 



Relativistic Fireball: Many workers have concluded that the relativistic fire- 
ball model is the most likely scenario for GRB creation and most prefer conver- 
sion of the blastwave energy into emission through blastwave interaction with 
an external medium. The basic outlines of the fireball model are probably best 



summarized by Waxman [ 119 1 : 



A compact tq ~ 10 7 cm source releases an energy E comparable to that 
observed in 7-rays, E ~ 10 51 ergs, over a time t < 100 s. The large 
energy density in the source results in an optically thick plasma that 
expands and accelerates to relativistic velocity. After an initial acceler- 
ation phase, the fireball energy is converted to proton kinetic energy. 
A cold shell of thickness ct is formed and continues to expand with 
time-independent Lorentz factor r ~ 300. The GRB is produced once 
the kinetic energy is dissipated at large radius, r > 10 13 cm, due to 
internal collisions within the ejecta and radiated as 7-rays through syn- 
chrotron and possibly inverse-Compton emission of blastwave accelerated 
electrons. Following internal collisions, that convert part of the energy 
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to radiation and that result from variations in r across the expanding 
shell, the fireball rapidly cools and continues to expand with approxi- 
mately uniform Lorentz factor r. As the cold shell expands, it drives a 
relativistic blastwave into the surrounding gas. 

Huang et al.|30| described this interaction at larger radius as: 

After producing the main GRB, the cooling fireball is expected to ex- 
pand as a thin shell into the interstellar medium (ISM) and generate a 
relativistic blastwave, although whether the expansion is highly radiative 



116 or adiabatic is still controversial. Afterglows at longer wavelengths 



are produced by the shocked ISM. 

Rees & Meszaros |8(| described the possibility of generation of GRBs from 
the interaction of a highly relativistic (r ~ 10 3 ) blastwave interacting with a 
tenuous ISM or CSM to generate a ~ 10 51 erg 7-ray burst. Meszaros et al.|fT0| 
worked this out in more detail (see also jnj for a review of fireball models). 
When the first optical afterglow of a GRB was finally found for GRB970228, 



although the data were very limited Wijers et al. [ 127 1 found support their model 
of an expanding, relativistic fireball of total energy ~ 10 51 ergs plowing into a 
surrounding medium. 

As early as 1993, Paczynski & Rhoads |8(| pointed out the similarity of GRB 
relativistic fireball interaction with CSM or ISM material to relativistic jets 
interacting with intergalactic material, supernova remnants (SNRs) interacting 
with the ISM, or RSNe interacting with the CSM and predicted the presence of 
delayed radio emission from the GRB process. Meszaros & Rees |72| and Vietri 
[ |116| also predicted late time radio, optical and soft X-ray emission based on the 
fireball model. 

Huang et al. |Q proposed that the extremely large r factors of the GRB- 
producing relativistic fireball decay quickly, within 3 to 4 days, so that the 
blastwave/CSM interaction generating the synchrotron emission is only slightly 
relativistic with 2 < r < 5. Even more rapidly evolving is the fireball model of 
Dai et al.j^ who proposed that the fireball is only mildly relativistic after ^3 
hours and non-relativistic after a few days. 



Waxman [ 1 19 1 studied the optical and radio afterglow from GRB970508 and 
found agreement with the cosmological fireball model for a fireball energy of 
~ 10 52 ergs and the blastwave expanding into an ambient medium with density ~ 



1 cm 3 . Waxman et al. [ 120 also compared the radio light curves of GRB970508 



to the predictions of the Waxman [ 1 1 9 1 fireball model and, although their fit is 



extremely poor after ~4 weeks of age, they ascribed the deviation to transition 
of the fireball from relativistic to subrelativistic expansion with a Lorentz factor 
r < 2 or, possibly, to a lower energy "jet" structure for the fireball at early 
epochs. Their model also predicts transition to a more uniformly emitting fireball 
structure by age ~5 weeks, becoming a spherical, non-relativistic emitter by age 
~12 weeks. 
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Relativistic Jets and Hypernovae: Support for the rclativistic fireball model 
is not universal, however. Dar p(| proposed an origin for GRBs and their longer 
wavelength afterglows in relativistic jets produced by mergers and/or accretion- 
induced collapse of compact stellar objects. He felt that such a scenario can better 
solve the "energy problem" (where some estimates of relativistic fireball energies 
can range up to a rather high > 10 52 ergs), the short time variability of many 
GRBs, and other apparent deviations from the relativistic fireball models. A 
somewhat different variant on the fireball model has been proposed by Paczynski 
[ jsi) who suggested that a "hypernova" explosion produces a "dirty" fireball from 
the violent death of a massive star. 



Multiple Origins: Livio & Waxman j6jj took the broader stance that the 
GRBs with afterglows, although probably not originating in neutron star-neutron 
star (ns-ns) mergers, could well have several possible origins including the col- 
lapse of massive stars as in SNe (GRB980519 and GRB980326) or could be 
produced by Black Hole-He star mergers (GRB990123 and GRB990510), with 
a third group such as GRB970228 and GRB970508 being possibly produced by 
Black Hole-He star mergers but without collimated jets being formed by the 
collapse. They disagreed with other workers that the transition from relativis- 
tic expansion to non-relativistic expansion of a spherical fireball will produce 
a detectable break in the afterglow decline rate. They also felt that the GRBs 
creating detectable afterglows are diverse with only some producing highly col- 
limated jets, but with all originating from massive stars and all interacting with 
a stellar wind environment. 

Perhaps the broadest summary and comparison of models with radio data for 
several GRBs previous to this review is that given by Chevalier & Li jl(| (see also, 
[ p"5| ) who were able to describe rather well the radio light curves of several GRB 
radio afterglows as arising from blastwave interaction with an external medium. 
They also suggested that there is evidence for two groupings: (a) GRBs arising 
from the explosion of a massive star, possibly a Wolf-Rayet star, embedded in a 
dense p oc r~ 2 stellar wind-formed CSM and accompanied by a supernova, and 
(b) GRBs arising from the explosion caused by a compact binary merger that 
occur in a constant density ISM and are not accompanied by a supernova. 



3 Summary 

We have assembled an overview of the radio emission from GRBs and have shown 
that there are many similarities to the much better understood radio emission 
from RSNe. In particular, we have shown that the formalism described by Weiler 



et al.[126| can be applied to the radio and OIR emission from a number of GRBs 
and provides a satisfactory description of their light curves over extremely broad 
frequency and time ranges. These results then imply that at least some GRBs 
probably arise in the explosions of massive stars in dense, highly structured 
circumstellar media. Such media are presumed to have been established by mass- 
loss from the preexplosion stellar systems. 
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